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Abstract. We have obtained the universal conductance distribution of two-dimensional disordered systems 
in the strongly localized limit. This distribution is directly related to the Tracy- Widom distribution, which 
has recently appeared in many different problems. We first map a forward scattering paths model into a 
problem of directed random polymers previously solved. We show numerically that the same distribution 
also applies to other forward scattering paths models and to the Anderson model. We show that most of 
the electric current follows a preferential percolation-type path. The particular form of the distribution 
depends on the type of leads used to measure the conductance. The application of a moderate magnetic 
field changes the average conductance and the size of fluctuations, but not the distribution when properly 
scaled. Although the presence of magnetic field changes the universality class, we show that the conductance 
distribution in the strongly localized limit is the same for both classes. 

PACS. 7 1.23.-k - 7 2.20.-i 



1 Introduction 

The hypothesis of single-parameter-scaling (SPS) [TlfS] con- 
stitutes the main foundation of our understanding of local- 
ization in disordered systems. The original formulation of 
the SPS hypothesis focused on the average conductance, 
but it was soon realized that the full distribution function 
of the conductance should be considered. Then, according 
to the SPS hypothesis, the conductance distribution func- 
tion should depend on a single parameter [5], for example, 
the mean conductance (g) or the mean of the logarithm. 

The validity of the SPS hypothesis has been thor- 
oughly checked in one-dimensional (ID) systems, where 
it has been shown that all the cumulants of In g scale lin- 
early with system size [T . Thus, the distribution function 
of In g approaches a Gaussian form for asymptotically long 
systems, which in ID are always strongly localized. It is 
fully characterized by two parameters, the mean (lug) and 
the variance of In 5. In the scaling regime, both param- 
eters are related to each other through a universal law. 



where L is the system size and ^ the localization length, 
defined as 

o T 

e = - hm — — . (2) 

L^oo (lug) 



Eq. was first derived in Ref. [5j within the so-called 
random phase hypothesis, which assumes that there exists 
a microscopic length scale over which phases of complex 
transmission and reflection coefficients become completely 
randomized. This relation reduces the two parameters of 
the distribution to only one and provides, therefore, a jus- 
tification and interpretation for SPS in ID systems. 

The situation in higher dimensional systems is not as 
clear as in ID systems. In those dimensions is far more 
difficult to do analytical calculations and numerical sim- 
ulations have been limited until recently to small sample 
sizes. In the diffusive regime, although the size of universal 
conductance ffuctuations depend on the dimension of the 
system, the distribution function of the conductance tends 
to a gaussian for all dimensions. Some people thought that 
in the localized regime it could happen something similar. 
In the strong localization regime, In g was claimed to be 
normally distributed in dimensions higher than one [6l[7l 
IHl. However, it was pointed out that the distribution is 
not log-normal [^ [TU1[TT1[T2] . We found numerically that 
the variance behaves as [T3|[T4] 

(T^ ^ A{-lng)" + B (3) 

with the exponent a equal to 2/3 in two-dimensional (2D) 
systems. In three-dimensional systems, preliminary results 
indicated a possible value of a = 2/5 ^14], but recent 
more extensive simulations point out to a value of 1/2, 
in agreement with directed polymer simulations [15]. The 
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constants A and B are model or geometry dependent. The 
precise knowledge of the dependence of cr^ with (— In g) 
made much easier the numerical verification of the SPS 
hypothesis, which we checked for the Anderson model [13^. 

In this paper, we concentrate in the strongly localized 
regime in 2D systems. Although experimental measure- 
ments of coherent transport at low temperatures are diffi- 
cult in this regime, knowledge of the conductance distribu- 
tion is of interest to better understand variable range hop- 
ping conductance, the metal-insulator transition in three 
dimensions or the crossover between the diffusive and the 
localized regime in 2D systems. 

In the strongly localized regime, the contribution of 
each Feynman path to the tunneling amplitude between 
two sites decays exponentially with its distance, and we 
can expect that some properties of the conductance (in 
particular the size dependence) should be dominated by 
the shortest or forward-scattering paths (FSP). This ap- 
proach was introduced by Nguyen, Spivak and Shklovskii 
(NSS) in a model to account for quantum interference 
effects in the localized regime. Medina and Kardar [17] 
studied in detail the model. They computed numerically 
the probability distribution for tunneling and found that 
it is approximately log-normal, with its variance increas- 
ing with distance as r^^^ for 2D systems. This is in con- 
trast with the ID case, where the variance grows linearly 
with distance, and with the implicit assumptions of some 
works on 2D systems. In our opinion, this approximation 
did not receive as much attention as it deserves, probably 
because in the SPS regime the localization length must be 
much larger than the lattice constant, which means that 
contributions from other paths cannot be negligible. 

A major step forward in the field of random systems 
was done by Tracy and Widom (TW) 18J who obtained 
the distribution function of the largest eigenvalue of ran- 
dom matrices belonging to the gaussian ensembles, or- 
thogonal, unitary and symplectic. It soon became clear 
that the TW distribution for the unitary ensemble also 
appears in the calculation of the length of the longest 
common subsequence in a random permutation |19j and 
in many other seemingly unrelated problems. For an in- 
troduction see Ref. [50]. Particularly important for our 
problem is the study by Johansson [21] of a specific type 
of directed polymer model. He obtained the distribution 
function of the lowest energy state exactly in terms of the 
Tracy-Widom (TW) distribution. Also relevant for us is 
the one-dimensional polynuclear growth model |22j , which 
studies the height fluctuations of a growing interface and 
is closely related to the Kardar-Parisi-Zhang equation [23 . 

We showed that for the 2D Anderson model in the 
strongly localized limit the conductance distribution is re- 
lated to the TW distribution [2^ . Here we present results 
for two different models in the FSP approximation, which 
is important to link exact results in directed polymers with 
the more realistic Anderson model. We show numerically 
that for a model with disorder with both positive and 
negative values (which cannot be exactly mapped to the 
solvable model) the skewness of the conductance distri- 
bution tends without free parameters the TW value. The 



role of percolation versus interference is also analyzed. We 
finally study the conductance distribution function in the 
presence of a magnetic field. We show that systems in 
the orthogonal gaussian ensemble (Anderson model with- 
out magnetic field) and in the unitary gaussian ensemble. 
(Anderson model with magnetic field) tend to the same 
behavior in the strongly localized limit. 

In the next section, we describe the two models that 
we have used in our calculations. In section III, we present 
a mapping of the conductance in localized systems to the 
free energy of directed polymers and use Johansson's re- 
sults to get the distribution function. In sections IV and 
V, we show that the distribution function of the logarithm 
of the conductance in the FSP approximation and in the 
Anderson model, respectively, are related to the TW dis- 
tribution. In section VI, we analyze the effects of a mag- 
netic field. We finalize with a discussion and conclusion 
section. 



2 Model 

We have studied numerically the Anderson model and 
its FSP approximation for 2D samples. For the Anderson 
model, we consider square samples of size Lx L described 
by the standard Anderson Hamiltonian 

-ff = ^ e^ajai -I- a]a,j + h.c. , (4) 

where the operator aj (a^) creates (destroys) an electron 
at site i of a cubic lattice and is the energy of this 
site chosen randomly between {—W/2, W/2) with uniform 
probability. The double sum runs over nearest neighbors. 
The hopping matrix element t is taken equal to —1, which 
set the energy scale, and the lattice constant equal to 1, 
setting the length scale. All calculations with the Ander- 
son model are done at an energy equal to 0.01, to avoid 
the center of the band. 

We have calculated the zero temperature conductance 
g from the Green functions. The conductance g is pro- 
portional to the transmission coefficient T between two 
semi-infinite leads attached at opposite sides of the sam- 
ple 



where the factor of 2 comes from the spin. ^From now on, 
we will measure the conductance in units of 2e^ /h. The 
transmission coefficient can be obtained from the Green 
function, which can be calculated propagating layer by 
layer with the recursive Green function method [251126] . 
This drastically reduces the computational effort. Instead 
of inverting an x matrix, we just have to invert L 
times L X L matrices. With the iterative method we can 
easily solve square samples with lateral dimension up to 
L = 400. We have considered ranges of disorder W equal 
to 13, 15 and 25, which correspond to localization lengths 
of 1.12, 2.4 and 3.2, respectively, and lateral dimensions up 
to L ~ 200 for the calculation of the distribution function. 
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Fig. 1. Schematic picture of the sample and the leads con- 
sidered: (a) narrow leads and (b) wide leads. The open cir- 
cles represent sites in the system and the solid dots sites in 
the leads. The lines represent the hopping between sites. The 
dashed lines are just a guide to the eye. 



which requires a huge number of independent runs to get 
good statistics in the tails. For this purpose, we average 
over a number of realizations larger than 6 x 10^ for each 
disorder and size. The leads serve to obtain the conductiv- 
ity from the transmission formula in a way well controlled 
theoretically and close to the experimental situation. 

To consider different possible geometries, we have used 
two types of leads: wide leads with the same width as 
the lateral dimension of the samples and narrow (one- 
dimensional) leads. These are attached to the sample at 
the centers of opposite edges, as shown in Fig[T|a). The 
scheme of the wide leads is shown in Fig[l|b). In both 
cases the leads are represented by the same hamiltonian 
as the system, Eq. Q, but without diagonal disorder. The 
narrow leads can be viewed as a simplified model of a point 
contact, while the wide leads should roughly correspond 
to electrodes in contact with the whole edge of the sample. 
We use cyclic periodic boundary conditions in the direc- 
tion perpendicular to the leads. 

The introduction of a uniform magnetic field B per- 
pendicular to the sample leads to complex hopping ma- 
trix elements. In the Landau gauge the vector potential 
is A = (0,— i?a;, 0). Then, the hopping matrix elements 
in the X direction are unchanged by the presence of the 
field, while the elements in the Y direction have to be mul- 
tiplied by the factor exp(±ia;_B), where the sign depends 



on whether we are connecting a site with the upper or 
lower site in the same column [57^. 

We have also studied the FSP approximation, first con- 
sidered by Nguyen, Spivak and Shklovskii |16j and widely 
studied by Medina and Kardar [17) in 2D square samples. 
One can write the matrix elements of the Green function 
between two sites a and b in terms of the locator expansion 



{a\G\b)=J2Y[ 



1 



r ier 



E 



(6) 



where the sum runs over all possible paths connecting the 
two sites a and h. In general the convergence of this series 
is very problematic, but in the strongly localized regime 
for distances much larger than the localization length one 
expect that the previous sum is dominated by the FSP. 
Considering only directed paths is well justified in the 
strongly localized regime, where the contribution of each 
trajectory is exponentially small in its length. We expect 
that back-scattering paths renormalize the site energies, 
but they should be irrelevant in the renormalization-group 
sense in the strongly localized regime. Based on this idea, 
the FSP approximation only considers directed paths. In 
this approximation, we will consider E = Q and two types 
of diagonal disorder: i) can only take two values W 
and —W, chosen at random with the same probability, 
which was the model originally considered in Ref . |16j ; it) 
ti = max[|a:|, l]sign(a;), where x is chosen randomly in the 
interval {—W/2,W/2) with uniform probability. The rea- 
son to substitute the small disorder energies, je^j < 1 by 
±1 is to take partially into account the effects of backward 
paths, which for these sites are important, avoiding at the 
same time problems of convergence in the locator expan- 
sion and ensuring that the transmission through any path 
is never larger than one. We will refer as NSS model to 
the FSP approximation with the first type of disorder. For 
the FSP approximation, we concentrate on the transmis- 
sion amplitude between two points in opposite corners of 
a square sample and we assume that the quantum trajec- 
tories joining these two points have to follow one of the 
(many) shortest possible paths. For the NSS model it is 
standard to use the path length /, rather than the system 
size L. For our geometry, I = 2L. The transmission at zero 
energy is equal to P7] 



T ■ 



21 



J\l) 



(7) 



where the transmission amplitude J{1) is given by the sum 
over all the directed paths 



directed 



J{1)= J2 Jr 



(8) 



The contribution of each path, Jp, is the product of the 
signs of the disorder along the path F. J does not depend 
on W for the first type of disorder and only weakly for the 
second type, for which we take W = 10. 

The variance of In 5 is entirely determined by J'^{1) and 
so it depends on /. It is convenient to quantify the mag- 
nitude of the fiuctuations in terms of the path length / in 
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the FSP approximation, (lug) is, of course, proportional 
to I, but the constant of proportionality depends on the 
disorder W. The assumption of directed paths facilitates 
the computational problem and makes feasible to handle 
system sizes much larger than with the Anderson hamil- 
tonian. The sum over the directed paths can be obtained 
propagating layer by layer the weight of all the trajectories 
in a very efficient way |17j . 



3 Mapping of the conductance in localized 
systems to polymers 

We have proven that for a specific FSP model the dis- 
tribution function of the logarithm of the conductance in 
2D systems is a Tracy- Widom distribution. To do so, we 
have to map the problem of the conductance in strongly 
localized systems to the problem of finding the energy of 
a polymer in a random environment. Each FPS path con- 
tributing to the conductance corresponds to a directed 
polymer. The calculation of the quantum amplitude be- 
tween two points, given by Eq. ([6]), in the FSP approxi- 
mation is then formally similar to the calculation of the 
partition function of directed polymers in a random po- 
tential at a very small temperature T = l/fc/3 

Z = Y,(^xp\-(3Y.h,\ = ,Y,l[exp{-ph,} (9) 



r ier 



hi are random site energies and F runs over all possi- 
ble configurations of the directed polymer. We can easily 
make Eqns. ^ and ([9| equivalent by associating (3hi with 
hi{E — ei). In this case, we can map the distribution of In 17 
in our system to the distribution of the free energy in di- 
rected polymers. As in polymer physics we consider real 
disorder energies hi, we then choose all the values E — ei 
to be positive. 

The distribution function of the ground state energy H 
of a polymer in a disordered environment was obtained ex- 
actly by Johansson [21 . In his model the random site en- 
ergies take integer values with probabilities Vvihi = k) = 
(1 — p)p'^. The ground state energy for polymers running 
between the origin and the point {x,y) is given by [2F 



H{x,y) 



2y/pxy + p{x + y) 

~P 
{pxy)^/^ 



(10) 



l-p 



{l+p) 




{x + y) 



2/3 



X2 



where X2 is a random variable with the TW distribution, 
corresponding to the distribution of the largest eigenvalue 
of a complex hermitian random matrix [18] . X2 verifies 



Pr(x2 <x)= f2{x')dx' = e-''^") 



(11) 



where g{x) is the solution of the equation 
g"{x)^u^{x), 



which tends to zero, g{x) —fQ,asx—f 00. u{x) is the 
global positive solution of the Painleve II equation 



= 2U-' 



(13) 



which tends to the Airy function, u{x) —f Ai(a;), when 
x 00. 

We can map Johansson's result for the polymer prob- 
lem to obtain the distribution function of the conductance 
for the localization problem by considering the FSP ap- 
proximation and that the disorder energies are of the form 



(14) 



with probability Pr(/c) = (1 ~ p)p^ for k = 0,1,2,- ••, 
in the limit /5 — s- 00. This is a very specific model, but 
we expect that their results apply in a much more gen- 
eral context, i.e., different disorder probabilities, as al- 
ready suggested by Johansson [51]. The same distribu- 
tion should also apply for non zero temperature, as in a 
renormalization-group sense temperature is an irrelevant 
parameter and the behavior is dominated by the zero- 
temperature fixed point. 

It is not trivial to extend the previous mapping to re- 
alistic problems with negative site energies, like the NSS 
and the Anderson model. It is then natural to study if 
the distribution function given by Eq. (11) also applies to 
these models. 



4 Conduction in the forward scattering paths 
approximation. 

First of all, we have studied the conductance distribution 
of the NSS model in order to check if it is also given by the 
TW function. For this model it is natural to consider In J^, 
which contains all the relevant information about fluctua- 
tions. In Fig.[2]we plot, for several system sizes, histograms 
of In as a function of z = (In — (In J^))/tT, where 
is the variance of In J^. The data are for narrow leads and 
each line corresponds to a different system size, specified 
in the figure. The thick solid line is the standardized TW 
distribution. We see how the results for the NSS model 
tend uniformly to the TW distribution as size increases, 
although the convergence is relatively slow. 

Although the results in Fig. 2 present noticeable finite 
size effects, it is clear that the TW distribution is fully con- 
sistent with the limiting distribution. This fact together 
with the previous mapping permit us to suggest that 



In J2 = + (3L^'^X2 



(15) 



in the limit L ^ 00. To better guarantee that the TW 
distribution is the proper limiting function, it is conve- 
nient to focus on the adimensional parameters of the dis- 
tribution, like the skewness (Sk = k^/ky^) or the kurtosis 
(Kur = ki/k"^). If the previous equation is correct, these 
two parameters must tend to the TW values as L — > 00 



independently of the values of a and (3. From Eq. ( 15 1 we 



(12) expect that the leading correction for the skewness, with 
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Fig. 2. Histograms of In as a function of the standardized 
variable z for several sizes of the NSS model. The thick solid 
line corresponds to the TW distribution. As size increases, the 
data approach the TW distribution. 

respect to the limiting value at L — > cxd, is proportional 
to L~^/^. In Fig. [3] we plot the skewness as a function 
of L^'^l'^ for the NSS model (circles) and for our second 
FSP model, with continuous disorder (triangles). The hor- 
izontal line corresponds to the TW value. We see how the 
skewness approaches the correct value and how it does so 
in the expected way for the two FSP models considered. 
Finite size effects for the NSS model are larger than for 
the other FSP model considered. A similar behavior is also 
observed for the kurtosis. It is clear from these trends that 
both models tend to the same universal distribution. 




0.00 0.05 0.10 0.15 0.20 



Fig. 3. Skewness vs L"^/^ for the NSS model (black circles) 
and the FSP approximation with continuous disorder (green 
triangles). The horizontal line corresponds to the TW value. 



We note that Johansson's results imply that the con- 
ductivity, in our case, is likely to be dominated by the 
most important FSP. On the contrary, the NSS model 
was designed to maximize the interference effects and all 



trajectories have the same amplitude. So, considering only 
the most conducting path has no sense in this case. It is 
necessary then to understand how the NSS may end up 
reproducing the TW distribution. In Fig. |4]we plot, for a 
given sample, the value of from the left corner to any 
other point in the sample. All points in a vertical slice 
are at the same hopping distance from the left corner. For 
each vertical slice the maximum value of is plotted in 
black and the minimum in white. Intermediate values are 
plotted on a gre y scale proportional to J^. This method 
gets rid of the exponential decay of with distance. How- 
ever, it does not take into account that the current must 
flow through the end point. Despite the fact that all single 
trajectories have exactly the same weight, interference ef- 
fects produce that most of the current is carried through 
very few well defined paths. The interference effects play a 
role only at short scales, producing a kind of renormaliza- 
tion process, which ends up in a dominant "renormalized" 
path. Probably, this renormalization also explains why fi- 
nite size effect are larger in the NSS model than in the 
other FSP model considered. 




Fig. 4. Plot of from the left corner to any other point 
in a given sample on a grey scale. On each vertical slice, black 
(white) corresponds to the maximum (minimum) conductance. 



5 Conductance for the Anderson model 

Once we have shown that the conductance in the FSP 
approximation follows the TW distribution, we study the 
Anderson model, which cannot be directly mapped to the 
polymer problem. We have calculated numerically the con- 
ductance for the 2D Anderson model and we have obtained 
its distribution function. By similarity with the FSP mod- 
els we first consider the narrow leads geometry. In Fig. [5] 
we plot histograms of In g for this model as a function of 
X = (In .9 — A)/i3, where A and B are chosen in order 
to have the same mean and variance as the theoretical 
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distribution, the TW distribution /2(x2)- The data are 
for several sizes and ranges of the disorder, W = 25 and 
L = 100 (squares), = 13 and L = 200 (dots), W ^ 15 
and L = 150 (triangles) and W = 9 and L = 70 (con- 
nected circles). The solid line corresponds to /2(X2)- All 
the cases in the strongly localized regime, represented by 
unconnected solid symbols, are fitted fairly well by /2(X2)- 
The agreement extends over more than four orders of mag- 
nitude. The connected empty circles correspond to a sys- 
tem near the crossover regime, whose distribution is close 
to a log- normal, represented by a dashed line in Fig. [5] For 
narrow leads, /2(X2) fits the data better than the gaussian 
function for L/^ > 13. 




X2 

Fig. 5. Histograms of Ing versus the scaled variable x for 
several sizes and disorders of the Anderson model with narrow 
leads. The continuous line corresponds to the TW distribution 
and the dashed line is a gaussian with the same mean and 
variance as the TW distribution. 



Considering the size dependence of the mean and the 
variance of In g and the excellent agreement between our 
data and the TW distribution, we conclude that in the 
strongly localized regime with narrow leads 



In g = — 



2L 
T 



7^ 



1/3 



X2 



(16) 



where /3 is a new constant and X2 a random variable with 
the TW distribution. In the SPS regime /? is a constant, 
independent of the disorder, the system size or the Fermi 
energy. We found, from present results and previous calcu- 
lations on the behavior of the variance, that it is approx- 
imately equal 3.4 for the Anderson model with narrow 
leads. Eq. (16) must be also valid outside the SPS regime 



(when the localization length is very small or when the 
Fermi energy lies in the band tails), but with a non uni- 
versal value of p. The data set in Fig. |5] for T4^ = 25 is 
outside the SPS regime, since it corresponds to a localiza- 
tion length of the order of the lattice spacing (^ = 0.97), 
and fits the TW distribution very well. The value of P is 
3.5 in this case. 



Table 1. Mean, variance, skewness, and kurtosis for the dis- 
tributions /2 and fo- 



/2 



fo 



Mean -1.77109 

Variance 0.81320 1.15039 

Skewness 0.2241 0.35941 

Kurtosis 0.09345 0.28916 



In table I we present the characteristic parameters of 
the distribution of /2(x2)- The fact that it has a mean 
different from zero, (X2) = —1.77109, has some practical 
consequences. According to Eq. ( 16 1, this implies a contri- 
bution to (In (7) proportional to L^^, already observed by 
us |12| in the Anderson model with narrow leads. This is 
of practical interest for the calculation of the localization 
length. Neglecting this term can cause errors of the order 
of 20 % in the estimate of the localization length. 

We did not find any L^^^ contribution for wide leads, 
which constitutes a strong indication that the conduc- 
tance distribution may depend on the leads, even in the 
strongly localized regime. The type of leads may change 
the universality class. We can expect a situation similar 
to polynuclear growth models ^22j, where the height dis- 
tribution was found to depend on the initial conditions. 
Our problem with narrow leads is directly related to the 
droplet model in Ref. [22_, which starts from an initial 
preferential point. With wide leads we have translational 
invariance and all the initial (and final) points are equiva- 
lent, a problem similar to stationary growth. In this case, 
the height fluctuations are described by the function fo 
whose accumulated distribution is [22; 

fo{x')dx' - [1 - (x + 2f" + 2g")g'] e-^^+^f^ (17) 

where f(x) is the solution of the equation 

f'[x) = -u{x), 
which tends to zero when x 



(18) 



00. g{x) and u{x) are 
given by Eqs. \12\ and (|13|, respectively. The character- 



istic parameters of the distribution /g are given in table 
I. As required by our previous finding that there were no 
L^/^ contributions to (In g), /q has zero mean. In Fig[6]we 
show the histograms of In g for the Anderson model with 
wide leads for several disorders and sizes, 14^ = 13 and 
L = 100 (squares), = 13 and i = 200 (dots), W = 15 
and L ~ 100 (triangles) and W — IQ and L = 50 (con- 
nected circles). The solid symbols correspond to systems 
in the strongly localized regime, while connected empty 
symbols to a case near the crossover. We plot in terms of 
X= (Iny — yl)/_B, where A and B are chosen in order to 
have the same mean and variance as the distribution /q. 
The full Une corresponds to /o and the dashed hne to a 
gaussian with the same mean and variance. The agreement 
between the numerical data and the theoretical distribu- 
tion is again excellent, showing that Ing satisfies in this 
case 

2L JL^^'^ 



Ing 



(19) 
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where /? is a new constant and xo ^ new random variable 
with distribution fa^x). In the SPS regime, /3 is a univer- 
sal constant approximately equal to 2.2. We note that, in 
the data for = 10 and L = 50, the right tail of the 
distribution is cut at xo ~ 4. This is due to the existence 
of a kink in the distribution for a transmission equal to 
one. For wide leads the distribution /o(xo) fits the data 
well when L/^ > 6. 




Fig. 6. Histograms of Ing versus the scaled variable x for 
several sizes and disorders of the Anderson model with wide 
leads. The continuous line corresponds to /o and the dashed 
line is a gaussian with the same mean and variance as fo- 

Our results support the idea that the Anderson model 
and its FSP approximation in 2D systems in the strongly 
localized regime will verify an equation of the form ( 16 1 



or ( 19 1 for any range of parameters, type of disorder, ge- 
ometry or boundary conditions. The distribution of the 
random variable will depend on boundary conditions, as 
we have shown for two particular cases, narrow and wide 
leads. It becomes clear that the definition of a "univer- 
sality class" requires to take into account the boundary 
conditions. A detailed study of these was done in the con- 
text of a PNG model [55]. By analogy with this model, we 
think that the concept of universality class remains valid 
as it seems that there exist only a finite set of universal 
distributions. It is difficult to know the complete set of 
universal distributions and to which of them will tend a 
general complex type of lead. 

Despite the existence of several universality classes, 
there are some results which are very robust, like the lo- 
calization length and the exponent of 1/3 in Eqs. ( 16 ) and 
( |T9| ) characterizing the size of fluctuations. This exponent 
determines the behavior of the cumulants and the tails of 
the distribution [28 . For the two types of leads studied 
we have 



lllfj{x) 



for 



oo 



(20) 



with do — 2 and d2 = 4. Eq. (20) is valid for the high 



five to boundary conditions, although both distributions 
/o and /2 behave as |22] 



In/, (a;) 



12 



for X 



'OO. 



(21) 



As we have mentioned, we consider L/^ > 6 as the 
condition for the strongly localized regime, where the con- 
ductance distribution is fitted by Eqs. ( 16 1 and ( 19 1 much 
better than by a log-normal. For mesoscopic samples close 
to the condition L/^ « 6 it should be possible to test ex- 
perimentally our predictions, since g ~ go exp(— 12) where 
go is close to one in units of 2e^//i. Our results can also 
be verified experimentally through the behaviour of the 
cumulants of the distribution. Eqs. (7) and (9) predict 
universal values for the skewness, kurtosis, etc, of the dis- 
tribution, given in table I. This limiting values can be ob- 
tained from a measurement of the cumulants in any range 
of parameter in the localized region, since Eq.(6) it is fairly 
well verified even near the crossover. From the tendency of 
the second and third cumulants, for example, it is possible 
to derive the asymptotic value for the skewness, which for 
wide leads should be ^3/^2^^ — 0.359 (see table I). 



6 Distribution function with a magnetic field 

We have studied the distribution function of lug in the 
localized regime of the Anderson model at zero tempera- 
ture when we apply a magnetic field perpendicular to the 
sample. The field changes the mean, producing negative 
magnetoresistance, and the variance of the distribution 
but not the distribution itself when it is rescale in terms 
of the variable x— (In <? — ^)/C, where A and C are again 
chosen in order to have the same mean and variance as the 
theoretical distributions /2 and /q. We found that the dis- 
tribution functions for both narrow and wide leads are the 



same as without a magnetic field, Eqs. (16 1 and (19 1 



In Fig. [7] we plot the distribution function of In g for 
two values of the disorder W = 15 and W = 25, two 
system sizes and several values of the magnetic field B 
(see the legend in the figure). We have considered the 
two types of leads used before: solid symbols correspond 
to wide leads and empty symbols to narrow leads. The 
solid lines are our theoretical distributions /2(X2) (nar- 
row leads, left curve) and /o(xo) (wide leads, right curve). 
We can see that all the points are fitted pretty well by 
the distributions /2 and fo, respectively, as in the absence 
of a magnetic field. We note that we have used periodic 
boundary conditions through out the paper, including this 
section. 

The results indicate that the magnetic field does not 
change the percolative nature of the conduction in the 
strongly localized regime. 



The value of /3 in Eqs. ( 16 1 and ( 19 1 changes very little 



conductance tail only. The other tail might be more sensi 



with the applied field, less than our uncertainty in the 
measurements . 

We have also checked that the application of a mag- 
netic field in the NSS model with narrow leads produces 
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Fig. 7. Histogram of Ingr for the Anderson model with an 
applied magnetic field for narrow (empty symbols) and wide 
(solid symbols) leads for the values of the disorder, system size 
and magnetic field given in the figure. The continuous lines 
correspond to ,f2(x2) and /o(xo)- 



negative magnetoresistance, but does not change the dis- 
tribution function. To be quantitative, we have plotted 
in Fig. 8 the skewness as a function of for several 

values of the magnetic field, detailed in the legend of the 
figure. The horizontal line corresponds to the skewness of 
the distribution /2 and the lower curve to the situation 
in the absence of a magnetic field. The skewness clearly 
tends to the TW value, as the system size tends to infin- 
ity, for all cases considered. We note that in the presence 
of a field the convergence is faster than in its absence. It 
is also interesting to note the dependence of the slopes in 
Fig. 8 with magnetic field. They increase with the field, 
while the value in the absence of a field is larger than in 
the presence of it. 



(/5 




Fig. 8. Skewness as a function of L"'^^'^ for several values of 
the magnetic field in the NSS model with narrow leads. The 
horizontal line corresponds to the TW value. 



7 Discussion and conclusions 

Present results confirm our previous belief that, in the 
strongly localized regime, directed path models are in the 
same universality class as the Anderson model [HllTSlfTT] . 
While the NSS model pretended to maximize interference 
effects, Johansson's model only considers the most impor- 
tant path. The agreement between both models indicates 
that it is percolation and not interference the dominant 
effect in this regime. We expect that the main effect of in- 
terference between different paths is a renormalization of 
the disorder energies. This information may be relevant to 
deal with interacting systems, since for many properties 
FSP models are a good approximation to the Anderson 
model and they can be extended to many-particle systems, 
whose direct simulation may be feasible for relatively large 
systems. 

The distribution functions of the conductance in the 
FPS approximation and the Anderson model appear in 
many other problems. Our results are fully consistent with 
the strong version of the SPS j29j if we take into account 
that boundary conditions (leads) are not irrelevant vari- 
ables in the renormalization group sense, and may change 
the universality class. We have checked that the skewness 
and the kurtosis obtained numerically for narrow leads 
tend to the theoretical predictions. We finally showed that 
the conductance distribution does not change when we ap- 
ply a magnetic field. 

The situation in 3D systems is more complex, since 
analytical approaches valid for 2D systems do not apply. 
There are no hints of distribution functions for any similar 
type of problem. Nevertheless, the quality of the fit in 2D 
systems suggests that numerical simulations will provide 
useful information. 
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